Human Genetics

Standard polygenic risk scores (PRSs) are constructed based on additive genome-wide association study (GWAS) summary statistics. Nonlinear machine learning methods have been increasingly applied to construct PRSs directly from individual-level data, with the aim of improving predictive performance over standard PRSs through their ability to model non-additive genetic effects. However, their superiority across studies has been inconsistent, and the conditions under which they provide meaningful improvements remain unclear. We combined theoretical analysis, simulations and a real-world application to investigate when two widely used nonlinear machine learning methods, random forest and XGBoost, outperform standard PRSs. Theoretical analysis showed that standard PRSs can implicitly capture part of the genetic variance attributable to nonadditive genetic effects through their contributions to marginal SNP effects, thereby losing less information than commonly assumed. Although nonlinear models have a higher theoretical potential, their greater flexibility incurs a bias-variance trade-off that can limit predictive gains at finite sample sizes. Simulations showed that XGBoost outperformed the standard PRS only when the genetic architecture involves a sufficiently large proportion of interaction genetic variance concentrated across relatively few interaction effects and large training samples were available. Random forest consistently underperformed the standard PRS. In an application to ischemic heart disease prediction using UK Biobank data, XGBoost showed no meaningful improvement in predictive performance over the standard PRS, whereas random forest again performed worse. Together, these findings suggest that nonlinear machine learning do not uniformly outperform standard PRSs; rather, their relative performance depends jointly on genetic architecture and training sample size. Our study helps to reconcile the inconsistent results reported across previous studies and provides a framework for identifying settings in which more complex PRS models are likely to be beneficial.

Sub-Saharan African populations carry pharmacogenomic alleles poorly represented in the European-derived reference panels underlying most clinical genotyping tools. We present a curated, machine-readable catalog of nine actionable alleles across six pharmacogenes (CYP2D6, CYP2B6, CYP2C9, CYP2C19, CYP3A5, NAT2) with African-specific frequency ranges, functional annotations, and evidence levels derived from reanalysis of 661 high-coverage whole-genome sequences across seven 1000 Genomes Project African populations. Direct comparison against PharmCAT v3.4.0 shows that CYP2D6 produces zero diplotype calls (0/661 samples callable) due to monomorphic reference positions absent from standard variant-only VCF output, a known limitation whose consequences for African allele carriers had not been reported. afripharmagen's reduced-position strategy identifies 243 CYP2D617 and 134 CYP2D629 carriers from the same input. For CYP2B6, CYP2C9, CYP2C19, and NAT2, both tools show concordance of 95-100%. Frequency gradients (CYP2B66: 30-50%; CYP2D617: 15-35% in West Africa; CYP3A5*1: 60-95%) translate directly into prescribing risk for efavirenz, tramadol, tacrolimus, and isoniazid. Pharmacogenomic decision support in African settings must incorporate population-specific allele definitions and input-format-aware strategies.

Obesity is a growing public health challenge, with body mass index (BMI) influenced by both genetic and environmental factors. While the role of direct genetic transmission is well established, evidence for genetic nurture effects, in which parental genotypes impact offspring through the environment, has remained mixed. This study investigates direct genetic transmission and genetic nurture effects on BMI across ages, using parent-offspring trios and pairs from the Dutch Lifelines cohort study (N = 18,897 offspring, aged 8 to 67 years). We leveraged the latest multi-ancestry BMI polygenic score (PGS) to construct transmitted (PGS-T) and non-transmitted (PGS-NT) polygenic scores, where PGS-NT consists of parental alleles not passed on to offspring and serves as a proxy for genetic nurture. Linear mixed models showed a large effect of PGS-T on offspring BMI (Beta = 0.416, p < 0.001), corresponding to a 1.85 kg/m2 increase per SD increase in PGS-T. PGS-NT had a small but significant effect (Beta = 0.026, p = 0.013), consistent with a genetic nurture effect accounting for approximately 6.6% of the effect of direct transmission. Parent-of-origin analyses showed that maternal PGS-NT effects were larger than paternal effects. PGS-T interactions with age indicated that direct transmission effects increased in childhood and stabilized in adulthood, while PGS-NT effects remained stable across age. Our findings suggest that direct genetic transmission is the dominant influence on BMI, while results are consistent with small genetic nurture effects that are driven by the maternal side.

Transmission-blocking vaccines (TBVs) are a promising strategy to reduce malaria transmission by targeting parasite stages within the mosquito. However, parasite genetic diversity may limit vaccine efficacy. We used next-generation amplicon deep sequencing to identify non-synonymous single nucleotide polymorphisms (SNPs) in Pfs25 from 184 Plasmodium falciparum isolates from Senegal, Tanzania, Ghana, and Burkina Faso. Prioritized SNPs were introduced into P. falciparum via CRISPR-Cas9. For the G116C variant, gametocyte development was evaluated by microscopy and qPCR, and mosquito infectivity was assessed by SMFAs. We identified 26 SNPs, including 24 novel variants. Functional assays showed that the Pfs25 G116C mutation did not affect gametocyte development or exflagellation. SMFA showed no significant differences in oocyst prevalence or intensity between mutant and WT parasites. These findings highlight the importance of integrating genetic surveillance with functional validation to guide the development of effective transmission blocking interventions

Background Lifestyle interventions are central to obesity prevention and management, yet interindividual variability in response remains incompletely understood. Here, we leveraged genetically defined, distinct obesity endotypes to examine lifestyle-body mass index (BMI) associations across biological pathways. Methods In the UK Biobank, we analyzed 305,713 participants with partitioned polygenic scores (pPSs) representing 10 obesity endotypes. We evaluated interactions between endotype-specific genetic susceptibility and physical activity, diet, sedentary behavior, and sleep on BMI using multivariable linear regression. Primary findings were externally evaluated in the All of Us Research Program using Fitbit-derived lifestyle measures. Results Favorable lifestyle behaviors were associated with lower BMI for all obesity endotypes, but the magnitude of these associations varied significantly across endotypes. Higher endotype-specific pPSs strengthened the benefits of physical activity (7 endotypes), healthy diet (3 endotypes), nonsedentary behavior (5 endotypes), and adequate sleep (7 endotypes) on BMI. Distinct endotypes demonstrated the greatest responsiveness to different lifestyle domains, with the metabolically unhealthy endotype showing the strongest interaction with physical activity, metabolically healthy endotype with sedentary behavior, hypothalamic dysregulation endotype with diet, and hypoinsulin 2 endotype with sleep, corresponding to differences in BMI of 0.22-0.49 kg/m2 between the highest and lowest pPS deciles. These interaction patterns were consistent in the All of Us cohort. Conclusions Obesity endotypes modify the association between lifestyle behaviors and BMI, demonstrating that responsiveness to lifestyle behaviors is heterogeneous and pathway dependent. These findings provide a framework for precision obesity prevention by identifying individuals who may derive greater benefit from specific lifestyle interventions.

Autosomal dominant gain-of-function mutations in platelet-derived growth factor receptor beta (PDGFRb) cause overgrowth of the skeleton and other connective tissue in Kosaki overgrowth syndrome. However, the target cell type and signaling pathways underlying PDGFRb-driven overgrowth are unknown. Normal postnatal growth is controlled by pituitary-secreted growth hormone (GH), which activates the STAT5 transcriptional factor to upregulate insulin-like growth factor 1 (IGF1). To investigate the role of the GH-STAT5-IGF1 pathway in PDGFRb-related overgrowth, we generated mice with a PDGFRb gain-of-function mutation in skeletal and fibroblast lineages, which resulted in STAT5 activation and gigantism. Conditional deletion of Stat5ab in connective tissue lineages rescued skeletal overgrowth and keloid-like fibrosis in the skin. Conditional deletion of GH receptor (Ghr) did not rescue overgrowth, indicating the physiological activator of STAT5 is not required for overgrowth. However, deletion of Igf1, the STAT5 target gene, and its receptor, Igf1r, in connective tissue, rescued the overgrowth phenotype. These findings demonstrate a GHR-independent STAT5-IGF1 signaling pathway in mutant connective tissue cells, which mediates PDGFRb-driven overgrowth in mice and potentially in humans with similar PDGFRB mutations.

Confounding is a central challenge in observational studies. Here, we propose a framework for identifying confounders of two non-causally related traits by employing cross-trait pleiotropy analysis to detect genetic loci that affect both traits and multi-trait colocalisation to identify molecular phenotypes mediating these effects. We apply this approach to the analysis of C-reactive protein (CRP) - a non-specific marker of inflammation - and 10 inflammation-related cancers. In UK Biobank, higher pre-diagnostic CRP levels are associated with increased risk of multiple cancers, but bidirectional Mendelian randomization provides little evidence for a causal relationship. Cross-trait genetic analyses identify 92 loci with shared CRP-cancer effects including those with established roles in cancer and 50 novel loci such as RSPO3 (breast cancer) and GCKR (colorectal cancer). Integration with proteomic and single-cell transcriptomic data identified putative molecular mediators at 24 loci including plasma TLR1 levels in breast cancer and CD4+ T cell IRF5 expression in kidney cancer. Notably, 15 candidate effector genes encode targets of approved or investigational medications, including IL6, PDE4D, and CASP8, indicating potential opportunities for their repurposing for cancer prevention. The proposed approach provides a generalisable framework for leveraging non-causal phenotypic relationships to yield insights into disease mechanisms and therapeutic targets for disease prevention.

Pathogenic variants in RUBCN, encoding the Run domain Beclin-1 interacting and cysteine-rich domain-containing protein (Rubicon) have been implicated in autosomal recessive spinocerebellar ataxia 15 (SCAR15). However, the molecular mechanisms underlying disease pathogenesis remain poorly understood. Here, we report 18 individuals from 15 unrelated families harbouring biallelic RUBCN variants, who present with an aggressive neurodevelopmental disorder variably characterized by seizures, developmental delay, intellectual disability and movement abnormalities that cause regression, progressive brain atrophy and neurodegenerative features. Through functional characterization, we demonstrate that a subset of disease-associated putative truncating variants disrupt autophagy regulation. In Caenorhabditis elegans models, loss-of-function RUBCN variants result in an increased autophagic flux and impaired neuronal function, recapitulating key features in humans. Correspondingly, cellular assays reveal that nonsense and frameshift RUBCN variants lead to defective autophagy inhibition, underscoring a crucial role for RUBCN as a key negative autophagy regulator. Molecular dynamics simulations rank the eleven missense variants by structural effect, with p.Arg813Trp alone altering the target protein at both the local and the regional level and lying within the RAB7A-binding module that the truncating alleles remove altogether. Our findings establish and expand the RUBCN-related disorders as a clinically and molecularly distinct subset of autophagy-related diseases. By delineating both the genetic landscape and cellular consequences of Rubicon dysfunction, this study enhances our understanding of autophagy-related neurodevelopmental disorders and provides a foundation for future therapeutic investigations.

Migraine is a leading cause of disability, yet preventive treatment remains largely empirical despite the availability of several mechanistically distinct therapies. Genetic data can clarify mechanisms and therapeutic hypotheses when association signals are integrated with molecular and clinical data. We meta-analyzed migraine GWAS data from 12 European ancestry cohorts (206,893 cases and 2,093,175 controls) and four African ancestry cohorts (22,115 cases and 178,626 controls). We identified 311 lead variants in European-ancestry analyses and 316 lead variants in trans-ancestry analysis. Fine-mapping and transcriptome-wide analyses prioritized variants and genes implicated in sensory neuronal signaling, vascular tone, and immune regulation, with convergent evidence at several established loci including TRPM8 and PHACTR1. Drug-repurposing analyses identified therapeutic targets and compounds, including established migraine treatments and candidates requiring experimental validation. Genetic correlations, Mendelian randomization, and a phenome-wide scan linked migraine liability to psychiatric, pain, and gastrointestinal phenotypes. Together, these findings expand the known genetic architecture of migraine across ancestries and provide a genetics-led map connecting association signals with biological pathways, multimorbidity and candidate therapeutic mechanisms, providing a foundation for future functional and translational studies.

Obesity is a significant public health concern. Early-onset obesity in the context of rare disease can reflect genetically-mediated pathology or elevated susceptibility through indirect mechanisms. Mapping the diverse characteristics and needs of young people with obesity in the rare disease population is a first step toward mechanistic and translational research. We carried out a retrospective comparative analysis of demographic, genotypic, phenotypic and health service utilisation data for young people with obesity (cases: n=500) and without obesity (controls: n=11,444) from the UK 100,000 Genomes Project rare disease cohort. Cases and controls were recruited prior to genomic diagnosis, across clinical disorder categories. We observed significant association between socioeconomic deprivation and obesity risk. Young people with obesity had significantly higher utilisations of acute care and mental health services, indicating an overall higher health burden. A curated panel of 519 candidate obesity-associated genes demonstrated aggregate association with obesity, although no single gene reached significance. Phenotypic comparison between cases and controls highlighted increased multi-organ and neurological system involvement, highlighting the overlap between neurodevelopmental and obesity risks. Within the case group, we conducted cluster analysis to identify early-onset obesity groups with different phenotypic profiles, potentially arising from different causal pathways - this identified six obesity subgroups of interest, with differing involvement of neurodevelopmental and other systems. Our study confirms that obesity co-occurs with a wide range of factors within the rare disease population, and is associated with significant physical and mental health needs, requiring holistic lifelong care.

Genotype imputation remains essential for large-scale human genetics studies, but its performance is limited by the size and ancestral diversity of available reference panels, reducing accuracy for rare variants and underrepresented populations. Here, we present a cloud-based imputation service built on a multi-ancestry reference panel derived from 515,579 jointly phased genomes from the All of Us (N=414,830) and National Human Genome Research Institute's Analysis, Visualization, and Informatics Lab-space (AnVIL, N=100,749) datasets. The All of Us + AnVIL reference panel is highly diverse and includes 261,163 participants most genetically similar to non-European reference populations, spanning 665,398,839 high-quality autosomal sites, representing a nearly 50% increase over TOPMed, the previous largest imputation service. Across multiple ancestry groups, the panel enables accurate imputation (empirical R2 0.8) for variants with allele frequencies as low as 0.2%, extending reliable imputation into the rare-variant frequency spectrum, including allele frequencies down to 0.002% and 0.006% for samples with European ancestry and African ancestry in the United States, respectively. Compared with TOPMed, the panel improves imputation accuracy across all ancestry groups except Africans, and recovers additional trait-associated variants not represented in existing reference panels. To facilitate broad community access while preserving participant privacy, we deploy the panel through a secure cloud-based imputation platform using privacy-preserving recombined haplotypes. This resource establishes a new foundation for genome-wide association studies (GWAS) and fine-mapping, especially in previously underrepresented populations.

Recent large-scale studies have enabled new knowledge about genetic underpinnings of morphological and electrophysiological alterations of the retina. Variation in retinal traits, often of neurodevelopmental origin, have been linked to major psychiatric disorders (MPDs). Here, we investigate the genetic overlap between MPDs and key retinal traits to identify underlying molecular mechanisms. We obtained genome-wide associations studies data for bipolar disorder (BD), major depression (MD), schizophrenia (SCZ), and the retinal traits retinal nerve fibre layer thickness (RNFL), ganglion cell inner plexiform layer thickness (GCIPL), and vertical cup-disc ratio (VCDR). We estimated the number of trait-influencing variants shared between traits with MiXeR and identified shared genetic loci with condFDR. Subsequently, we examined the biological pathways of the genes mapped to shared loci. This revealed that GCIPL shared the most genetic variants with MPDs (~60%), followed by RNFL (~40%), and VCDR (~20%). The genetic variants shared between retinal traits and MPDs showed disorder-specific patterns with more pronounced overlaps of SCZ and BD with RNFL, and MD negatively correlated with GCIPL. Gene-pathway analysis highlighted the importance of GABAergic neurotransmission and a two-stage neurodevelopmental process in SCZ, whereas the role of mitochondria and a weaker developmental component were observed in BD. The results also implicated synaptic functioning and gene-expression processes in MD. Furthermore, polygenic analysis suggested that the genetic architecture of retinal traits can distinguish between MPDs. Our findings indicate shared genetic underpinnings between retinal traits and SCZ, BD, and MD, implicating altered neurodevelopment and neurotransmission underlying the retinal link to major psychiatric disorders.

Unmeasured confounding threatens causal inference and replicability in observational multi-omic studies across variable environments. Genetic instrumental variables (Mendelian randomization) and negative-control calibration each address complementary sources of unmeasured confounding, yet no existing framework unifies them for omics-scale mediation analysis. We introduce ICONIC, an R package that embeds genetic instruments and negative controls within a proximal causal inference framework for total-effect and mediation analysis. ICONIC implements eight estimators spanning five confounding-control strategies, supports continuous, binary, and time-to-event outcomes, and provides extensive diagnostics including sensitivity analyses that map estimator performance across plausible assumptions. Ground-truth benchmarks are calibrated to real-omics covariance structures via a hybrid generative model (GAN + feature-level Gaussian copula) rather than parametric simulation, and a companion planning tool predicts performance gains from collecting additional omic data. We demonstrate ICONIC in two case studies: identifying placental transcriptomic mediators of gestational diabetes on birth weight (n = 164), and tumor-expression mediators of smoking intensity on lung cancer survival (n = 494). Notably, ICONIC's diagnostics recommended different estimation strategies across the two scenarios, reflecting differences in the likely influence of unmeasured confounding. ICONIC is freely available at https://github.com/sbresnahan/iconic/.

Chemoradiotherapy (CRT) is the standard-of-care therapy for many solid malignancies, yet predictive biomarkers of treatment response remain limited. We identified a germline single nucleotide polymorphism (SNP) in an intrinsically disordered region of the lysine demethylase KDM3C/JMJD1C (p.S464T) that is associated with CRT outcomes in locally advanced rectal cancers (LARC) and head and neck squamous cell carcinoma (LA-HNSCC). In silico modeling with AlphaFold predicted S464T substitution influenced interaction between phosphorylated KDM3C and RNF8 FHA domain. In cellular models, conversion of S464 to T464 increased sensitivity to DNA-damaging agents. S464T substitution impaired damage-induced MDC1-RAP80 signaling and downstream RAP80-BRCA1 colocalization. SNP carrying cells impaired DNA repair causing genotoxic stress that is associated with increased cGAS-cGAMP innate immune signaling and increased apoptosis. Population analyses with the SNP highlighted an increase incidence of UV-induced skin and other cancers, linking inherited variation in the chromatin regulatory gene KDM3C to genome instability, cancer risk, and therapeutic vulnerability.

In primary healthcare centers lacking advanced imaging, community-based active tuberculosis (TB) case finding often relies on basic symptom screening. This approach often misses cases and leads to the inefficient allocation of rapid molecular testing (RMT). We aimed to develop and internally validate a simple clinical triage scorecard to improve TB detection and guide RMT use in resource-constrained settings. We conducted a retrospective cross-sectional study of 15,137 adults ([&ge;]18 years) evaluated within the Zero TB Yogyakarta program (2020-2025). Participants with complete clinical assessments and confirmatory GeneXpert results were included. Using multivariable logistic regression, we identified independent clinical predictors, which were subsequently transformed into an integer-based point scorecard. Model performance was evaluated via discrimination and calibration, utilizing bootstrap resampling (1,000 iterations) for internal validation. Among the 15,137 participants, 251 (1.7%) were GeneXpert-positive. The final multivariable model identified eight independent predictors: age, male sex, body mass index, prolonged cough, hemoptysis, unexplained weight loss, TB contact history, and diabetes mellitus. The model demonstrated strong predictive accuracy, with an optimism-adjusted AUROC of 0.836 and good calibration. When translated to the integer scorecard and compared directly to standard national symptom screening, the scorecard performed (AUROC 0.81 vs. 0.73; p<0.001). At a high sensitivity cut off score of [&ge;] 0, the tool achieved 93.63% sensitivity and 41.33% specificity. This point-of-care clinical scorecard provides higher diagnostic accuracy than standard symptom screening algorithms. By offering flexible operational thresholds, it empowers local health programs to dynamically balance the urgency of case detection with available diagnostic capacity, optimizing GeneXpert allocation where advanced radiological imaging is unavailable.

Background: Many at-risk patients lack access to genetic services due to a genetic counselor (GC) workforce shortage. Little is known about how digital alternatives impact patients with and without cancer who meet criteria for genetic testing. Methods: eREACH2 is a randomized 4-arm non-inferiority trial where pre-test (visit 1) and/or return of results (visit 2) GC counseling was replaced with a patient-centered digital intervention. Arms include: A (GC/GC), B (GC/digital), C (digital/GC) and D (digital/digital). Primary outcomes were non-inferiority in uptake of genetic services and change in genetic knowledge and general anxiety from baseline to post-disclosure of results (T0-T2). Secondary cognitive and affective outcomes were assessed using non-inferiority ANOVAs and equivalency chi-squared tests in intention-to-treat and per-protocol analyses. Findings: 773 participants were recruited nationwide; 46.6% from rural areas. Mean age was 51 years (range 20-87), 13% male, 12% non-white, 29% had less than a college education, and 33% had a personal history of cancer. 584 (76%) patients completed testing (14% had a positive result, 16% had a VUS). In the primary ITT analyses, we met the non-inferiority for uptake of genetic services and anxiety, but results were inconclusive for knowledge. Secondary outcomes were heterogeneous across arms. Arm C demonstrated consistently favorable effects, while Arms B and D showed less favorable outcomes in select domains (e.g. satisfaction and MICRA). Patients who received positive or VUS results via digital disclosure had significantly higher MICRA scores - indicating greater negative response to testing. Interpretation: In this large, randomized trial of patients with and without cancer, the eREACH intervention was effective for pre-test counseling, but inconclusive for digital disclosure of results. Exploratory analyses suggest that digital delivery could be a reasonable alternative for individuals receiving negative results, while those receiving positive or VUS results may derive some short-term psychosocial benefit from GC disclosure.

Purpose: To evaluate whether initiation of GLP-1 receptor agonists (GLP-1RAs) is associated with anti-VEGF treatment burden in type 2 diabetes patients with diabetic macular edema (DME) in the IRIS(R) Registry (Intelligent Research in Sight). Methods: Incident GLP-1RA initiators were matched 1:1 with controls via Mahalanobis distance matching (9,896 pairs; N=19,792) on sociodemographics, DME risk factors, and factors influencing GLP-1RA prescription including hypertension, obesity, chronic kidney disease. A longitudinal mixed-effects event-study model evaluated monthly anti-VEGF injection frequency over a 36-month window (12 months before through 24 months after initiation), adjusting for DME duration. Visual acuity (VA) and central subfield thickness (CST) were secondary outcomes. Results: Following GLP-1RA initiation, anti-VEGF injection trajectories did not significantly differ between the matched GLP-1RA and control cohorts (interaction coefficients -0.18 to 1.59, P>0.05). Likewise, no differences in VA were observed between cohorts (-0.05 to 0.04 logMAR, P>0.05) or CST (-14.12 to 33.58 {micro}m, P>0.05). Conclusion: In these matched cohorts, GLP-1RA initiation was not associated with the trajectory of anti-VEGF use or changes in VA or CST. Precis We used the American Academy of Ophthalmology IRIS(R) Registry (Intelligent Research in Sight) to identify patients with DME. In 19,792 matched patients, there was no significant reduction in injection frequency post GLP1-RA initiation and no significant change in VA or CST.

Background: Refractory and macrolide-resistant Mycoplasma pneumoniae pneumonia (MPP) has emerged as a major challenge in pediatric respiratory medicine, amplified by the post-2023 resurgence. However, a systematic overview of the research landscape specific to treatment-refractory and drugresistant disease in children remains lacking. Methods: Research articles and reviews on pediatric refractory or macrolide-resistant MPP published between 2000 and 2025 were retrieved from OpenAlex using Boolean searches. After screening, 2,286 records were quantitatively analyzed for annual output, contributing countries/institutions, thematic clusters, and citation-burst dynamics using Python. Results: Annual publications grew exponentially, with a pronounced surge after 2023 (n=378 in 2025). China produced the highest volume (45.1%) but recorded fewer citations per publication than the US, Japan, and Canada. The literature resolved into four clusters: macrolide resistance/molecular basis, epidemiology, etiology/co-infection, and refractory disease management. Burst analysis showed an evolution from earlier fronts like 23S rRNA mutations and azithromycin to recent emerging trends like pandemic-related co-circulation, genotype surveillance, and co-infection. Conclusions: Research on pediatric refractory and resistant MPP is expanding rapidly, shifting in emphasis from etiologic descriptions toward resistance mechanisms and clinical management. Standardizing the treatment of macrolide-unresponsive disease and post-pandemic epidemiological surveillance represent the principal directions for future work. Keywords: Mycoplasma pneumoniae; children; macrolide resistance; refractory pneumonia; bibliometric analysis; research trends

Motivation: Cell-free DNA methylation provides a minimally invasive signal for early cancer detection and tissue-of-origin prediction. Most methods represent methylation measurements as independent fixed-window features and therefore do not explicitly model relationships among genomic regions. Results: We developed PANGEM (Pan-cancer Graph-based Cancer Detection Using the Cell-free DNA Methylome), a graph-learning framework that represents genomic bins as nodes and integrates CpG context, genomic proximity, and sample-specific methylation similarity in the graph topology. Across five repeated stratified train-test splits, PANGEM achieved the highest mean performance among evaluated methods, with an AUROC/AUPR of 0.997/1.000 for binary cancer detection and macro-AUROC/AUPR of 0.977/0.870 for multiclass tissue-of-origin prediction. In the independent INSPIRE cohort, 72 of 78 cancer cases (92.3%) exceeded the binary classification threshold, and PANGEM correctly classified 9 of 17 head and neck cancer cases (52.9%), the highest accuracy among evaluated methods. Subnetwork analysis further identified recurrent, graph-connected methylation patterns, including a 111-DMR subnetwork with increased methylation in cancer samples.

Single-cell technologies enable linking disease-risk variants to gene regulatory effects in specific cell-state contexts. However, most so called "single-cell eQTL" studies use a "pseudobulking" strategy to identify expression Quantitative Trait Loci (eQTLs), obscuring subtle dynamic regulatory effects of disease alleles. Here, we propose Dynema (Dynamic eQTL mapping in single cells) for fast and accurate genome-wide mapping of context-dependent and independent eQTL effects at true single-cell resolution. To identify eQTLs, Dynema uses a Poisson model with cluster robust variance estimators (CRVEs) to account for correlation of single-cell profiles from the same individual. In contrast to other common methods, Dynema achieves statistical calibration and scales to genome-wide analysis in large single-cell datasets in realistic timeframes. We applied Dynema to two independent T cell datasets and identified reproducible cell-state-dependent eQTL effects. Some cell-state-dependent eQTLs are missed by pseudobulking approaches, and many others are conditionally independent from lead eQTL effects. We show that TSPAN32 and other autoimmune loci colocalize with cell-state-dependent eQTLs. Mapping context-dependent eQTLs at single-cell resolution enables the definition of the molecular effects of complex disease alleles.